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Abstract 



We study the Cauchy problem for the Korteweg-de Vries (KdV) hierarchy in the 
small dispersion limit where e — >■ 0. For negative analytic initial data with a single 
negative hump, we prove that for small times, the solution is approximated by the 
solution to the hyperbolic transport equation which corresponds to e = 0. Near 
the time of gradient catastrophe for the transport equation, we show that the solu- 
tion to the KdV hierarchy is approximated by a particular Painleve transcendent. 
O ^1 This supports Dubrovins universality conjecture concerning the critical behavior of 

Hamiltonian perturbations of hyperbolic equations. We use the Riemann-Hilbert 
approach to prove our results. 



1 Introduction 



^ , In this manuscript, we will prove a particular case of a conjecture in [T3] about the 

CNj I formation of dispersive shocks [26j in a class of Hamiltonian perturbations of the quasi- 

O ■ linear transport equation 

^ ' ut + a{u)ux = 0, 

where u = u{x, t), x,u £ M, t G and where a is an arbitrary regular function of u. 
We will restrict ourselves to the case of the KdV hierarchy, which, for any m G N, is a 
Hamiltonian perturbation of the equation 

><■ Ut^ + CmU'^Ux = C^ = (-l)"'+^ ^"'^^"',+ ^^" . (1.1) 

The equations in the KdV hierarchy can be written in the form 

ut^-{-irdxiJm{u,eux,e\xx,...,e^'^d^x"'u) = 0, G R+, m G N, (1.2) 

where e > and dx = The function ipm is polynomial in its variables and it is the 
variational derivative of the Hamiltonian Tim \20\ l22l I40j 

ipmiu,eux, ■■■) = f^". , ^-171= [ hm{u,euxix),e^uxx{x), . . .)dx, 

which is defined as 

SV-m _ dhmiu,€Ux,. . .) _ d dhmiu,€Ux, ■ ■ ■) ^ d^ dhm{u,eUx, ■ ■ .) _ 
6u{x) du dx dux dx"^ duxx 

The Hamiltonians 1-Lra satisfy the Lenard-Magri recurrence relation [33] 

Slim n^^m-l . , ^ r , [ U^{x) 

^^M^ = ^ M^' - = j (1.3) 

P = e^dl + Audx + 2ux. 
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The densities are uniquely determined by together with the conditions 

u t 2 \| \ ^) '^rn. m+2 
hra{u,eu^,e Ua;x,---)U=0 = 1 — TTT -TTV^ • 

(m + l)(m + 2) 

For m = 1 the equation (|1.2p coincides with the standard KdV equation 

ut^ + 6uux + e'^Uccxx = 0, e > 0, (1.4) 
and for m = 2, 3 one has the equations 

Ut2 — SOu^Ux — {20UxUxx + WuUxxx) — ^^Uxxxxx = 0, (1-5) 

+ UOu^Ux + (70t(^ + 280uUxUxx + TOn^ns^) 

+ e"^ {70u2xU3x + ^2uxU4x + 14iiii5x) + e^ujx = 0. (1.6) 

The corresponding Hamiltonians for m = 1 , 2 are |2H [23} [36] 

Each equation in the hierarchy can also be written in the Lax form |30] 

Lt^ = [L,Am], (1.7) 

where L is the Schrodinger operator 

L = e'^dl + u, (1.8) 

Lf^ is the operator of multiplication by ut^ , and Am is an antisymmetric higher order 
operator with leading order (_i)™+i4™'£2m^2m+i^ The lower order terms are deter- 
mined by the requirement that [L, Am] is an operator of multiplication with a function 
depending on u, Ux, ■ ■ ■ , U{2m+i)x- For m = 1 and m = 2, we have 

Ai = 4:e'^dl + 3{udx + dxu), (1.9) 
A2 = -I6e^dl - 20e\d^iu + uO'^) + 5{dxUxx + Uxxdx) - l^{dxU^ + u^dx). (1.10) 

We will study the behavior of solutions to the KdV hierarchy in the small dispersion 
limit where e — )• 0. 

When e = the KdV hierarchy reduces to the transport equation (|l.ip . Let us 
assume that we have sufficiently smooth negative initial data uq[x) with a single local 
minimum, and which tend to rapidly at ±00. The Cauchy problem for (jl.ip can then 
be solved implicitly using the method of characteristics, which leads to 

u{x,tm) = Uo{0, -X + CmUo{irtm + i = ^- (I-H) 

This describes a left-moving solution for any m G N, but the part of the solution near 
the minimum moves faster than the less negative parts, so that the slope steepens at 
the left of the negative hump as tm increases. The time tm where the slope becomes 
vertical is called the time of gradient catastrophe, and is given by 

C = ^- zi • (1-12) 

mmax^eKl-C^-w™ (0^o(0} 

After this time, the solution to (jl.lip is no longer single- valued. 
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In order to determine the point of gradient catastrophe x'^ and u'^ = u{x'^,t^), one 
has to solve the system of three equations 

F{u; X, t^) := -x + CmU^'tm + hiu) = 0, (1.13) 
F'{u; X, tm) = mC„^u"'-H^ + /^(n) = 0, (1.14) 
F"{u; X, tm) = m{m - l)C„n'"-2t„ + /'[{u) = 0, (1.15) 

for the unknowns u,x,tm- Here primes denote derivatives with respect to u, and 
is the inverse of the decreasing part of the initial data uq. Among the possibly many 
solutions of (|1.13p - (|1.15|) . the point of gradient catastrophe is the solution {u'^;x'^,tm) 
with minimal time t^- We say that the gradient catastrophe is generic if 

k := -F'"(n^;x^^^) = -m(m - l)(m - 2)C7„(n^)— - /["(u^) ^ 0. (1.16) 

In the case m = 1, it is well-known that the dispersive term e'^Uxxx in (|1.4p reg- 
ularizes the gradient catastrophe that occurs for the Hopf equation ut-^ + 6uux = 0: 
the solution to the KdV equation exists for all ti > under suitable conditions on the 
initial data uo{x). For ti < if , the KdV solution is approximated by the Hopf solution 
for small e > 0, and for ti > tf, an interval of rapid oscillations is formed where the 
KdV solution can be modeled using Jacobi elliptic 6'-functions [26l [3T1 [39| [TT | [T2 t [24]. 

For m > 1 we have not been able to find results about global existence in time in the 
literature. For initial data in weighted Sobolev spaces only local results stating that the 
solution exists (and stays in the same space) for small times tm seem to be available, 
see e.g. |28j . For the class of analytic initial data we will consider, global existence in 
time should not be an issue for any m E N. We will comment on this later, see Remark 
12.11 However our aim is not to prove an existence result for the KdV hierarchy, rather 
we want to study the asymptotic behavior for small e of the solution u{x,tm,^) to 
the Cauchy problem for (jl.2p under the assumption that it exists. Before gradient 
catastrophe, one does not expect to see a significant difference between solutions to 
(jl.ip and (jl.2p for small e. After the gradient catastrophe, the KdV hierarchy solution 
will also develop a region of oscillatory behavior for small e. 

The purpose of this manuscript is twofold. First, for tm < tm, we will prove that 
u{x, tm; e) = 'u{x, tm) + O(e^) as e — >• 0, where u{x, tm) is the solution to the dispersion- 
less equation. Secondly, for tm ~ x ^ x'^, we will show that the KdV hierarchy solu- 
tion can be approximated for small e by a special Painleve transcendent U = U{X,T), 
which solves the fourth order ODE [2?! [2] fT3l l29] 



X = TU 



+ ^iU^x + 2U Uxx) + ^Uxxxx 



(1.17) 



This ODE is the second member of the Painleve I hierarchy, and we refer to it as the 
Pj equation. The relevant solution is real and has the asymptotic behavior 

U{X,T) = t{6\X\)^/^T-6'^/^T\X\-^/^ + 0{\X\-^), asX^ioo, (1.18) 

3 

for any fixed T E M, and has no poles for real values of X and T [SJ [Ml I35j. It is also 
remarkable that U {X, T) is an exact solution to the KdV equation normalized as 

Ut + UUx + ^Uxxx = 0. (1.19) 
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It was conjectured by Dubrovin in [13j that, for any Hamiltonian perturbation of a 
hyperbolic equation [l3l[32], a generic solution u{x,t,e) has an asymptotic expansion 
of the form 

u{x, t, e)=u^ + aie^/^f/ (a2e-6/^(x - x'' - a3{t - t^)), a4e-^/^(t - t")) + o (e^/^) , 

(1.20) 

for x,t near the point of gradient catastrophe x'^jt'^ of the unperturbed equation, and 
with constants ai, a2, 03, 04 depending only on the initial data and on the equation. The 
expansion should hold in a double scaling limit where e — )• 0, but at the same time x and 
t should tend to the point x'^ and time of gradient catastrophe for the unperturbed 
equation in such a way that the arguments of U in (jl.20p remain bounded. In other 
words, the limit is such that e — )• and at the same time e~^/^(x — x^ — 03(1 — t^)) and 
— t^) remain bounded. The Painleve transcendent U{X,T) is thus conjectured 
to describe the behavior of the solution to the perturbed equation near the point of 
gradient catastrophe for the unperturbed equation, and is expected to be universal 
in the sense that it is independent of the choice of the equation and independent of 
the choice of initial data. The only quantities in (|1.20p that depend on the initial 
data and on the equation are the constants 01,02, 03, 04, and the values of x'^,t^,u'^. 
The asymptotic formula (jl.20p was shown numerically for a certain class of equations 
including the KdV equation and the second member of the KdV hierarchy [25^ I15j. 

We prove this conjecture in the special case of the KdV hierarchy for a class of 
analytic initial data with a single negative hump. So far, the conjecture had been 
proven only for the KdV equation [3]. Similar results appear also in double scaling 
limits for Hermitian random matrix ensembles [6] and in the semiclassical limit of the 
focusing nonlinear Schrodinger equation |14j . 

1.1 Statement of results 

We study the Cauchy problem for equation (jl.2p with m G N. Similarly as in [3j, we 
impose the following conditions on the initial data uq. 

Assumptions 1.1 

(a) uq{x) is real analytic and has an analytic continuation to the complex plane in 
the domain 

S = {z gC : \lmz\ < tan6'|Rez|} U {z G C : |Imz| < a} 
where < 6 < 11/2 and a > 0; 

(b) uq(x) decays as \x\ — )■ 00 m 5 such that 

uo{x) = O (^^) , s > 0, X G 5, (1.21) 

(c) for real x, uq{x) < and uq has a single local minimum at a certain point xm, 
with 

Uq{xm) = 0, Uq{xa.i) > 0. 

Without loss of generality, we assume that uq is normalized such that uq{xm) = 
-1. 
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We prove the following result. 



Theorem 1.2 Let uo{x) satisfy the conditions described in Assumptions [Ll\ and let 
m G N, X G M and tm < f^, where is the time of gradient catastrophe for given 
by Ill.lS^} . If u{x,tm-,€) solves equation with initial condition n(x,0, e) = uq{x), 

then we have 

u{x,tm,e) = u{x,tm) + Oie"^), ase—^0, (1.22) 

where u{x, tm) is the solution to the Cauchy problem for M.l]) . given by \1.11\) . 

Remark 1.3 We will prove this result in detail for values of x,tm where u{x,tm) is 
decreasing, i.e. for x < xm + Cmtm since xm + Cmtm is the position of the minimum at 
time tm- For x > xm + Cmtm-, the proof is similar, but several technical modifications 
are needed. We will discuss those changes in Remark 14. II 

Our second result describes the behavior of u{x, tm, e) when x, tm approach the point 
and time of gradient catastrophe at appropriate speeds. We will prove that (11.20p holds, 
with the values of the constants given by 

as = Cm(y!^T, a4 = • (1.24) 

We will give an asymptotic expansion as e — t- for (x,tm) in a space-time window of 
size C'(e^/'^), and in addition x — Xc must be equal to Cm{u'^)"^{t — tm) plus a correction 
of size ©(e^/"^). In this shrinking (as e — t- 0) region in the (x, tm)-plane, the transition 
takes place between asymptotics determined by (jl.ip and the oscillatory asymptotics 
that are expected to be present for tm > f^. The transition is described by the 
Painleve transcendent U{X,T). In addition, we will also compute the next term in the 
asymptotic expansion (|1.20p . which is of order 0{e^^'^). This term is rather complicated 
but can still be expressed completely in terms of the Painleve transcendent U{X,T). 
We have an expansion of the form 

n(x, tm, e)=Uc + aie^l^U + cie^/^ {QUx + Uxx + - Scaf/r) 

+ Cie^''{2UxQT + ^UUt + hjxxT) + 0(e'/^), (1-25) 
where we used the abbreviations 

U = U fa2e-6/^(x -x"- asit - O), a4e-^/^(t - t^)) , (1.26) 



Q = ^.UxUxxx - % + - + ^ + ^^^i- (1-2T) 
The values of ci , C2 , C3 are 

mCm{u'^)"'-\tm-t'f,,) ( 2(m-l) 2F^A^\ 

= 4k^^ l^^^ + ^T^J' (^-^"^ 

with k given by ()1.16p . and F by (jl.lSp . 
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Theorem 1.4 Let uq{x) satisfy the conditions described in Assumptions [7T7] and as- 
sume that the generic condition 

k := -F"'{u^) = -m{m - l)(m - 2)K)"-3 - fl'{u^) / (1.31) 

holds with m G N. Write u'^,x'^,t^ for the solution to the system ^.13]) - ^1.14^ - fLT5\} . 
Let us take a double scaling limit where e — )■ and at the same time x — t- and 
im ~^ tm such a way that, for some X, T € M, 

hm e = X, lim ^ = T. (1.32) 

If u{x,tm,€^) solves equation hi. 2(1 with initial condition n(x,0, e) = uq(x), the asymp- 
totic expansion lil.25\) holds in the double scaling limit. 

The proofs of our results are based on the direct and inverse scattering transform 
for the KdV hierarchy. This approach rehes on the Lax representation (jl.7p . and the 
inverse scattering transform can be formulated as a Riemann-Hilbert problem, where 
one searches for a function which satisfies a prescribed jump condition, depending on 
the reflection coefficient for the Schrodinger equation Lf = A/. Solving the RH problem 
asymptotically as e — )• leads to small dispersion asymptotics for u{x, tm,e). We will use 
a Deift/Zhou steepest descent method similar to the one in [HI HI] for the asymptotic 
analysis of the RH problem. This method consists of a series of transformations M i— )• 
r I— 7- S" I— >• i? of the RH problem, which results at the end in a RH problem for R which 
can be solved approximately for small e. The first transformation M ^ T involves the 
construction of a G-function satisfying convenient jump and asymptotic conditions. 
The second transformation T ^ S deforms the jump contour from the real line to 
a lens-shaped contour. The last transformation S ^ R requires the construction of 
local and global parametrices. A local Airy parametrix will be needed for the proof of 
Theorem 11.21 For the proof of Theorem ll.4| we will need to build a local parametrix out 
of a model RH problem related to the Pj equation. The most important new features 
here compared to [3] are the generalization of the G-function to the case m > 1, 
and the generalized construction of the local Pj parametrix. 

2 Riemann-Hilbert problem for the KdV hierarchy 

We construct a RH problem using particular solutions to the Schrodinger equation 
Lf = A/, with L given by (jl.Sp with potential u = u{x, tm, e). This construction is well 
understood |101 [11 [38l El [16] , but we summarize the main lines here for the convenience 
of the reader. 

For negative u, the Schrodinger operator has no point spectrum. Moreover, if u 
solves the equation (jl.2p . as a consequence of the Lax equation (jl.7p . the eigenvalues 
of L are independent of t^. Since our initial data uq are negative, it follows that L has 
no point spectrum at any time tm > 0. If 

\u{x,tm, + x'^)dx < OO, (2-1) 

there exist [TO] fundamental Jost solutions V'± = V'i i^] x, tm, e) and 4)± = 4>± (A; x, tm, e) 
to the Schrodinger equation satisfying the asymptotic conditions 



lim ,^±(z;x,t^,e)^tv^" = 1, (2.3) 



lim ^±(z;x,t„,e)e±?^-^" = 1, (2.2) 
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for A € C \ {0}. We fix V— A to be the principal brancli of the square root which is 
analytic for A G C \ M"*" and positive for A < 0. We consider ^lJ± and (j)± as functions in 
the variable A, whereas x, tm, e will be parameters. Those solutions can be constructed 
as a solution to Volterra integral equations as in [10], and the analysis of the integral 
equations shows that tlj_ and can be continued analytically for A in the lower half 
plane, and and to the upper half plane. Also asymptotics as A — t- oo for (j)± and 
ip± can be deduced. 

The fundamental solutions ^p± and (j)± are related as follows, 

with 

|a(A)|2 - |6(A)p = 1, for A < 0, (2.5) 
and where a and b are independent of x. The quantities 



a{X;tm,€)' ' ' ' a(A;tm,e)' 

are the reflection and transmission coefficients (from the left) for the Schrodinger equa- 
tion and depend on tm,e through u{x,tm,€)- They are continuous for A < 0. In 
particular we have 

t±^[f^i^^^^ x^-oo, ^^^^^ 
a y^(i ^ ) X +00. 

If u solves the higher order KdV equation (|1.2p . the Gardner-Greene-Kruskal-Miura 
|22j relations are 

da db 2i ™ , , . 2m+i . , 

= 0, — = -4-t„(-A)^6. (2.7) 

Indeed, p.7|) implies that ^'4'+ + ^mV'+ is also a solution to the Schrodinger equation. 
The asymptotics as x — t- -|-c« then imply that 

— V'+ = -^^V'+ + (-A) 2 

(Xlm. £ 



Together with (j2.6p this implies (j2.7p and thus 



r(A;fm,e) = r(A;0,e)e 



„. 2m+l 



The transmission coefficient is analytic for A in the upper half plane, and the reflection 
coefficient is analytic in a region of the form {tt — 6*0 < arg A < tt}, with 6*0 > 0. If the 
potential u is smooth in x, the reflection coefficient decays rapidly as A — )• — oo. 

We will now construct the solution to a RH problem using the Jost solutions 'i\}± 
and write M = M{X;x,tm,f) by 



M{\]X,tr. 




(2i 
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where 



4>± = 4>±{^\x,tm,e), (j)± = cl)±{X]x,tm,e), (2.9) 



a = o(A;e), o*=a(A;e), (2.10) 

and 

a(A;x,t„) =:E(-A)5+4"^t^(-A)™+5. (2.11) 



Using (12. 4p and the asymptotics for the Jost solutions as A — )■ oo, one shows that M 
solves a RH problem, see [H [381 

RH problem for M 

(a) M(A;x,tm;e) is analytic for A G C\M, 

(b) M has continuous boundary conditions M±{X) for A G M \ {0} that satisfy the 
jump conditions 

^+(^) = ^-(^)U(A;.)e— \'-|.o(A;.)P J ^ < 

(2.12) 

M+(A) = M_(A)f7i, fTi=(^J forA>0, 

(2.13) 

with ro(A; e) := r(A; 0, e). 

(c) We have 



M{X;x,tm,e)= \ | ( I - —^^3 + ©(A'^) ) , for A ^ 00, 



(2.14) 



1 1 

^i\/—X —i\f—X/ 



with g = q{x, tm, e) independent of A, and u = is the solution to equation (jl.2p 
with initial data uq. 

In other words, the solution to the KdV hierarchy with initial data uq can be 
recovered from the solution to the RH problem by the formula 

d 

u{x,tm,e) = -2iedxMi^ii{x,tm,e), ^ i'^-'^^) 

where Mn{X; e) = 1 + ^•^i.ii(;Mm, e) ^ ^^^^^^ as A ^ 00. 



Remark 2.1 Using the vanishing lemma approach developed in [191 [THl [T7] , one can 
show that the RH problem for M is solvable for any value of x, f G M, e > if rg 
has sufficient regularity and sufficient decay at —00 and if |r'o(A)| < 1 for A < 0. 
The solvability of the RH problem can be used to prove that the Cauchy problem for 
equation ()1.2p is solvable for initial data in a suitable space following the proofs in [41] . 
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3 Asymptotic analysis of the RH problem as e — > 



This section contains the asymptotic analysis of the RH problem for M as e — )• 0. 
We will study the RH problem for fixed x,tm (for the proof of Theorem II. 2p . and 
in a double scaling limit for x ^ x'^, tm ^ (for the proof of Theorem ll.4p . The 
construction of the G-function in Section 13.11 is almost the same in both cases, and 
also for the transformations M i— )• T i— )• S" in Sections 13.21 and 13.31 and the construction 
of the outside parametrix in Section 13. 4| we do not need to distinguish between the 
regular and the critical case. It is only when we construct a local parametrix that there 
is an essential difference. For the regular case, we assume that x and tm are such that 
Uxix,tm) < 0, i.e. we consider only the decreasing part where x < xm + Cmtm, the 
position of the minimum at time tm- For the increasing part, several changes have to 
be made, see Remark 14.11 



3.1 The G- function 

In this section, we will define a G-function, which will be needed to modify the jumps 
of the RH problem in a suitable way. Let us first write the Abel transform 

r ^^ . ^ ^ F{(;x,tm)d^ x ^ r i m 

FA{X;x,tm) = 2 — — ' 

where F is defined by (jl.l3p . We define the G-function G = G(A; x, tm] u) as 



nf\ + \ Vu- X /° FA{v;x,tm) , 

G{X;x,tm;u) = / ^ dr]. (3.1) 

We will have to make two different choices for u G (—1,0). For the proof of Theorem 
11.21 we need to choose u = u{x,tm) to be the solution to equation (jl.ip : this solution 
is given by (jl.lip or, equivalently, by (|1.13p . For the proof of Theorem 11.41 on the 
other hand, we fix n to be n = n'^ = u{x^,t'f^). Whenever the choice of x, and tm is 
unimportant below, we will simply write Gm{X) for Gm{X;x,tm]u). 

G is analytic for A G C \ [n, +oo). As A — )• oo, we have G(A) = ©(A"^/^), and 
writing 

Gi{x,tm]u) := lim (-A)1/2g(A; x, t^, tx), (3.2) 

A— >oo 

one verifies that, for both choices of u made above, we have the identity 

u 

dxGi{x,tm) = -■ (3.3) 

For A on the branch cut, G satisfies the jump properties 

G+(A) + G_(A) = 0, for A G (0,+oo), (3.4) 

G+(A) + G_(A)-2/9(A) + 2a(A) =0, forAG(u,0), (3.5) 

with p given by 

The function p(X) is clearly well-defined for —1 < A < 0, but we can extend it to 
an analytic function. Because of the analyticity of the initial data uq, p is analytic 
in a neighborhood of (—1,0). It cannot be extended to an analytic function in a full 
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neighborhood of —1 because the inverse //, behaves hke a square root near —1, but we 
can extend it to an analytic function in a neighborhood of [—1,0), except for a branch 
cut which we choose along A < —1. We write = 17+ U U (—1, 0) for such a region 
where p is analytic, with 17+ C C+ and $7^ = 17+. For later convenience, we choose 
17+ sufficiently small so that it is contained in the region {tt — 9q < argA < vr} where 
ro(A) is analytic (see Fig. [1]) 




Figure 1: The regions 17+ and 17 . 

Let us define an auxiliary function (j) in 17+ U 17^ by 

(p{X; X, tm) = G{\; x, tm) - p{X) + a(A; x, t^), (3.7) 

so that 4) is analytic across (— l,u), but not on (—1 — (5, —1), because this is a part of 
the branch cut for p, and not on (ti, 0), because this is part of the branch cut for G. It 
is an analytic function for A G 17 \ [w, 0). By (j3.5p . we have 

2(^+{\;x,tm) = G+{X;x,t„,)-G-[\;x,tm), forAG(u,0). (3.8) 

After a straightforward integral calculation as in [3j , one observes that (j) can be written 
as 

^{\;x,t^) = (3.9) 

2 3 

= -Vu- XF{u;x,tm) + 3('"~ ^)^F'{u;x,tm) 

-^{u - X)If"{u; x,tm) ^'^^^^ " (^-^^^ 



3.2 First transformation M ^ T 

We are now ready to perform a first transformation of the RH problem. This will lead 
to more convenient jump matrices for the asymptotic analysis as e — )> 0. A crucial role 
is played by the G-function and its properties discussed before. 

Define 

r(A; X, t^, e) = {\ M(A; x, i^, e)e-7G(A;x,t™)<x3. (g.n) 
We then have 

RH problem for T 

(a) T is analytic in C \ M, 
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(b) r+(A) = r_(A)t>T(A), as A G M, with 
VTiX) = et^-W'^3^;A^(A)e-t^+W'^^ 
and VM is the jump matrix in the RH problem for M given by 



(3.12) 



-ro(A; e)e-2^"(Wm)A i _ \r^(^x. e)p 



, for A < 0, 



for A > 0. 



(3.13) 



(c) As A — )• oo, 



r(A) = (/ + o(A-i)) f. ^ 



iv— A — w— A 



(3.14) 



The left multiplication with the triangular matrix in ()3.1ip was needed to transform 
(|2.14p in (|3.14p . but has no effect on the jumps for T. The jump matrix (|3.12p can be 
simplified using the properties ()3.4p - ()3.5p of the G-function and the definition (13. 7p of 
(p. We write the jump matrix in a different form depending on the value of A G R. For 
A > 0, by (f3l4l) we have 

VTiX) = ai. (3.15) 

For a sufficiently small choice of 5i > 0, we can write the jump matrix in terms of the 
previously defined function (p on the interval (—1 — (5i,0). For A S {u,0), we use (j3.5p 
and (13.81) to conclude that 



ut(A) 



iK+(A) 



^ iK*_{\) (l-|ro(A)|2)eT<^+Wy 
where we have written k for 

K{X;e) = -iro{X;e)eTpW, for A G 17+, 



as A G {u, 0), 



(3.16) 



(3.17) 



with boundary values on M denoted by k_|_(A) and k!^(A) = k+(A). For A G (— 1 — (5i,n), 
by dHS]) and 



vt{X) 



1 m+(A)eT'/'+W^ 
MX)e-T't'-W l-|ro(A)|2 ^ 



(3.18) 



Here, the boundary values (j)± are needed only on (—1 — 6i, —1], on (— l,?i), we have 
(p = (p±. Finally, on (— oo, —1 — 6i), we have 



vt{X) 

Using 
equation by 



1 



ro(A)ef(e(A)+-(A))^ 



as A < -1 - 5i. (3.19) 



-fo(A)e-f (e{A)+"(A)) 1 _ |^^(;^)|2 J ' 

, (|3.1ip . and (|2.15p . we recover the solution of the higher order KdV 



u{x,tm,€) = 2dxGi{x,tm]u) -2iedxTl^{x,tm,() 
= u-2iedxTii{x,tm,e), 



(3.20) 
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where T^^ is given by 



Tu{X;x,tm,e) = 1 + ^'^^^/^''^ + 0{X-'), as A 



oo. 

A ■ 

The aim of this RH analysis is to end up with jump matrices that decay to the 
identity matrix when e — )■ 0. To get a feeling for the small e behavior of the jump 
matrix vt, we need to have information about the reflection coefficient ro(A;e) for 
small values of e. We have the following results, see Oil] in combination with [37J, for 
any choice oi Si > 0. 

(i) For A < -1-^1, 

ro(A; e) = ©(e-'^^/^), as e ^ 0, (3.21) 
with c > 0. 

(ii) for A lying in a region as defined in Section 13. H but A bounded away from 
— 1, say |A + 1| > ^, we have 

K{X;e) = l + 0{e), as e ^ 0, (3.22) 

(iii) for A G (n + (5i, 0), we have 

(1 - |ro(A)|2)eT'^+W = O(e-^A), c> 0, as e ^ 0. (3.23) 

The latter was shown in [3] for m = 1 and u = only, but the same argument 
applies to the case m > 1 and u = u{x,tm)- 

This implies that the jump matrix ^^(A) tends to iai for A G (w, 0), and that it 
tends to / for A < —1 — 5i. On (—1 — 5i,u), the jump matrix is oscillatory for small 
e. In the next section, we will deform the contour in such a way that the oscillatory 
behavior turns into exponential decay. 

3.3 Opening of the lens T S 

The jump matrix vt{X) can be written in the following factorized form for —1 — 6i < 

X < u, 

Because the first factor is analytic in a complex region il+ and the second in il.^, this 
factorization can be used to deform the jump contour: the interval (—1 — Si,u) can be 
deformed to a lens-shaped contour as shown in Figure [21 

Define S as follows, 



S{X) 




-iK{X)e 
1 



in region I, 



, , / 1 0\ (3.25) 

r(A), elsewhere, 

with K*(A) = k{X). Since k (resp. k*) is analytic in (resp. 1^+), S is analytic in 
each of the regions in Figure [2] if we choose the lens sufficiently close to the real line so 
that it lies in Q. 

These are the RH conditions for S. 
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Vt 




Figure 2: The jump contour S5 and the jumps for S 



RH problem for S 

(a) S is analytic in C \ Ss, 

(b) S+{X) = S^{X)vs for A G Ss, with 



on S 



^ fR(A) (l-|ro(A)|2)eT'/'+wl ' 



on E2, 



wt(A), 



as A € (ti, 0), 
as A G (-00, -1 - 61) U (0, +00). 



(c) S{X) = {I + 0{X~')) (^.J- 



-X -W-x 

For large A, S'(A) = T(A), and by (fM]) we have 
u{x, t,e) = u — 2iedxSli{x, t, e), 

where 

Sli{x, t, e 



as A 



00. 



Sii{X;x,t,e) = 1 + 



-A 



+ 0{X as A — )• 00. 



(3.26) 



(3.27) 



(3.28) 



For u = u{x,tm) in the definition of the G-function ()3.ip . we have F{u;x,tm) = if x 
belongs to the decreasing part of u{x,tm)- For tm < we also have F'{^;x,tm) < 
for ,f G (— li 0). If the contours Si and S2 are chosen sufficiently close to the real line, 
one uses ()3.9p to verify that, for any neighborhood U of u, there is a c > such that 



Im0(A; X, tm] u) > c, 
llll(j){X; X , tm] u) < — c, 
Im0+(A;x,tm;u) < -c. 



for A G Si \^/, 
A G ^2\U, 
for A G (u,0) 



(3.29) 
(3.30) 
(3.31) 



For u = = u{x^,tm), we have F{u]x^,tm) = and F'{S^]x^,tm) < for G (—1,0) \ 
{u}, and one again uses (j3.10p to check that the inequalities (j3.29p - (l3.3ip hold. Since 
(j) is continuous in x and tm, there must be a (5 > such that (j3.29p - ()3.3ip hold also 
for u = n^, and for |x — x'^l < 5, \tm — tm\ < S (possibly for a smaller c > 0). The 
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inequalities imply together with (j3.2ip . (|3.22|) . and (j3.23p that the jump matrix vs for 
S is uniformly close to constant matrices on \ as e — )• 0: we have exponential 
decay to / for A < —1 — 5i and on (Si U S2) \ U, and decay to iai on {u, 0) \ U, with 
an error of order 0(e). 



3.4 Outside parametrix 

We now deal with the jumps that do not tend to the identity matrix as e — )• 0. Therefore 
we ignore for a moment all jumps that tend to I as e — )• 0, and a small neighborhood 
of u. We are then left with a jump iai on (n, 0), and a jump ai on (0,+oo). We 
can explicitly construct a function with those jump properties, and with the same 
asymptotic behavior as S. Indeed, if we define 

p(-)(A) = (-A)i/^(n-A)-'^3/4 j^l (3.32) 
we have 

RH problem for p(°°) 

(a) : C \ [n, +00) ^ C2x2 ig analytic, 

(b) satisfies the following jump conditions on (u,+oo), 

pjoo) ^ p^oo)^^^ ^ ^ ^^^^ (3 



(c) p(°°) has the asymptotic behavior 

1 



/2 -i(-A)l/2 



, as A —7- 00. (3.35) 



This outside parametrix will determine the leading order asymptotic behavior of S as 
e — )■ 0, but to prove this, we need to control the jump matrices also in the vicinity of 
u. We will do this by constructing a local parametrix near u. 



3.5 Local parametrix near u 

So far, there was no need to distinguish between the regular case, where x, tm are fixed 
and where u = u{x,tm), and the singular case, where x,tm are involved in a double 
scaling limit together with e and where u = u'^ = u{x'^,t^). The choices of u will 
however be crucial for the construction of the local parametrix. In either case, we want 
to construct a local parametrix in a neighborhood Z// of u in such a way that 

RH problem for P 

(a) P : Z7 \ ^ C2x2 is analytic, 

(b) P satisfies the following jump condition on n S5, 

P+(A) = P_(A)^;p(A), (3.36) 
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Figure 3: The jump contour for P 



with vp given by 



U 1 



vp{X) 



<t>{X;x,tm) ^ j 

i 0. 



^ 2i 



as A G Si, 



1 0\ 

2i^M.„ ^ ^ \ 1 as A G 2.2, 



, as A G (n^, 0), 



(3.37) 



(c) As e — )• 0, we have the matching condition 

P(A) = (/ + 0(eT))p(°°)(A), forAeSW, 7 > 0, 
between the local parametrix and the outside parametrix. 



(3.38) 



Such a parametrix has, for e — 0, approximately the same jumps as S has onU r\ S5. 
The idea is that P will approximate S near u for small e, whereas will be a good 
approximation elsewhere. 

Note that by (j3.10p . (f>{u;x,tm) behaves like c{u — A)^/^ as A — )• n for tm < f^, 
but (l){u'^,x^,t^) behaves like c(u'^ — A)''/^ as A — u'^, so the jump matrix vp behaves 
differently in the regular case than in the singular case in the vicinity of u. This is the 
reason why we need to build different local parametrices in both cases. In the regular 
case, we will be able to construct a local parametrix with the value of 7 in ()3.38p equal 
to 1, in the singular case we will have a weaker matching with 7 = 1/7. 

3.5.1 The regular case: tm < tm 



Here we let u = u{x,tm) be defined by the equation (|l.lip . i.e. u{x,tm) is the solution 
to the unperturbed equation (II. ID . Then we will construct the local parametrix near u 
using the following model RH problem. 



RH problem for 

(a) ^ is analytic for C £ C \ F, with F 

271 

3 



271 1 

" 3 ) /■ 



{C G C : argC = 0,argC 



f-,orargC 
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(b) ^ satisfies the jump relations 



c^+(C) = <i>-(C) 
ci>+(C) = f-(C) 
$+(C) = <i>-(C) 



e-2e(c) i 
i 

1 ie2e(f) 
1 



for arg C = 0, 



(3.39) 



for arg C = ^. (3.40) 



1 



for aigC 



27r 
' 3 ' 



(3.41) 



where 6{() = IC^^^i with branch cut on (— cxd,0). 



(c) $ has the following behavior as C — ^ oo, 



-1 



i + o{c 



(3.42) 



with the branch cuts of (— C) " along (0, +oo). 



Several RH problems equivalent to this one have appeared in the literature, see e.g. 
[SI IHl El [H] , and it is well-known that this problem can be solved explicitly in terms of 
the Airy function. The unique solution to this precise problem has been constructed 
in H]. In the singular case later on, we will use the same model RH problem, but with 
a different value for 9. We only need the existence of a 2 x 2 matrix- valued function 
<I> satisfying the above RH conditions, the precise construction in terms of the Airy 
function is irrelevant. 



We define the parametrix P = P{X; x, tr, 
PiX) = E{X;eMe-'/^fiX)), 



by 



(3.43) 



where E is an analytic function in U and / is a conformal map from U to a neighborhood 
of 0. Let us first define / by the requirement that 



9{f{X)) = -f{Xf/' = ^{X), 



(3.44) 



and 



(3.45) 

By dHIO]), 
0, which is true because of 



(piX) = ±i(/)(A), for ±Im A > 0, 

so that the branch cut of (j) coincides with the one of /^^^ along X < u. 
equation (I3.44p only defines / analytically if F{u;x,tm) 
our choice of u. Then we have 

/(n) = 0, /'(n) = F'{u; x, 1^?/^ > 0. (3.46) 

Now we can choose the lens in U in such a way that / maps CiU to the jump 
contour T for Then P has jumps on Eg n Z// which are given by (I3.37p . 

Now define E by 



i?(A;e) = i=p(-)(A)(| "/)(-e-^/^ 



/(A))- 



(3.47) 



Then it is easily verified that E is analytic in U, and using ()3.42p . we obtain the 
matching 

P(A)P(°^)(A)-^ = P^^\X) {I + 0{e)) P(~)(A)-^ = / + 0{e), (3.48) 

for A € dU as e — )• 0. When x,tm approach x'^,t'^, the uniform convergence breaks 
down because f'{u) approaches zero. Therefore we will now construct a different local 
parametrix in the singular case. 
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3.5.2 The singular case: double scaling limit e — )• 0, x — )■ x^, tm — ^ 

In this case we need a slightly modified model RH problem, which is obtained by 
replacing 9 in the RH problem for ^> by 

r (C; X, T) = ^C'/2 - |c3/2 + XQ^'\ (3.49) 

We then obtain the RH problem 

RH problem for <1>'^ 

(a) is analytic for C e C \ with = {C G C : argC = 0,argC = ^,orargC = 
7 1 /• 

(b) satisfies the jump relations 

n(C) = ^-(C)(^ . qJ, forargC = 0, (3.50) 

^(0 = ^1(0(^0 \ J, forargC = f. (3.51) 

$;(C) = ^i(C)(.^2r(c) °), forargC = -f. (3.52) 

(c) has the following behavior as C — t- oo, 

*'(c) ^ ^(-c)-j" (_\ ;) (/+ .«.3(-c)-"^ - i («; (-0- 

In [5] , existence of a solution for a slightly different but equivalent RH problem has 
been proved, and the transformation to the RH problem for has been given in [S]. 
We write = $'^(C; T) for the solution to this RH problem. In the regular case, the 
constants Q, U vanished, but this is no longer true in the singular case: they now depend 
in a transcendental way on the parameters X, T. We have Q = Q{X, T), U = U{X, T), 
where U solves the Pj equation with asymptotics given by (jl.lSp . Furthermore we have 



the relations 

Qx = U, (3.54) 

Wx = Q'^U + U'^, (3.55) 

V = -QU-^Ux. (3.56) 



Those identities follow from the fact that <I>3s:*^'^'~ is a polynomial in C. Substituting 
(|3.53|) gives a function with terms proportional to C"^; . . . as C ^ co. The require- 
ment that those terms vanish, leads to (|3.54p - ()3.56p . Using the equation ()1.17p and the 
asymptotics for Q obtained in [5|, it is straighforward to derive the identity 
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For any Xq,Tq € M, there are complex neighborhoods of Xq and Tq such that the 
RH problem is solvable for X, T in those neighborhoods, and such that ()3.42p holds 
uniformly in X, T. 

We want to construct a parametrix of the form 

P(A) = ii;(A;e)cI>^(e-2/V(A);e-6/V(A;x,t„),e-^/^<72(A;t„)) . (3.58) 



where f,gi,g2,E are analytic functions which are to be determined. If we want P to 
satisfy the jump conditions given in (j3.36p . we need to construct /, gi, 52 in such a way 
that 

0" (e-2/7/(A);e-6/V(A;x,t„),e-4/752(A;t„.)) = ^0(A), (3.59) 
or equivalently 

(/(A); 5i(A; x, t^),g2{\; t^)) = ?(A). (3.60) 

Moreover we want that e^^^^5i(A; x, tm), e~^'''^52(A; fm) remain bounded in the double 
scaling limit where e — )• 0, x — )• j;'^, — ^ for fixed e > 0. This is the reason why we 
need to restrict to the scalings (jl.32p . 

If we take u = = u{x'^^ f^), we have 

F(n^;x^4) = F'(n^;x^4) = F"(u'=;x^4) = 0, 
and (j3.10p implies 

^(A; X, e) = -VX-W^iFiu"; x, t^) - F{u^; x^ 4)) 

-^(A-u^)t(r(n^;x,^^)-F'(u^;x^4))-^(A-n^)t(F"(u^;x,^^)-r'(n^;x^4)) 



-^J^^F"'{e,x,tm){X-0-'d^. (3.61) 



Then it is clear that the first line in ()3.6ip vanishes like a square root as A — )• u'^, 

the second line vanishes like c(A — u*^)^/^, and the third line behaves like c'(A — u^y^'^. 
Therefore we define /, gi , g2 by the equations 

lUs^^^^'^'^"^/'^'"^^'"^''*^^^^"^)'^^ (3.62) 

5i(A; X, tm)f{\f'^ = -^/A^(F(u^; x, tm) - F{u''- x\ 4)), (3.63) 

- ^^(^/(A)3/2 = -^(A-n^)i(F'(n-x,t„) -F'(n-x^4)) 
_ 1(A - u-)\{F"(u-- X, t^) - F"{u'- x\ 4)) 
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4 
15 



/^F'"(e;x,^^)-F'"(^;x^4))(A-0§^ie• (3-64) 
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One verifies that this defines /, 171 , §2 analytically in U and moreover in such a way that 
()3.60p holds. A straightforward calculation yields 

f{u^) = 0, /'K) = (8A;)7, /"(n^) = (3.65) 

gi{u'';x,tm) = 9'i{u'',x,tm) = -^^^^^gi{u'',x,tm), (3.66) 

g2[u ,tm) , (3.67) 

, ^ , /2(m-l) 3f"{u-) 



g'2{u';t^) = [-^^^ - JJJ^ ) 92{u'';tm), (3.68) 
where 

X = X - Cmiu'^rt^, k = -F"'{u^)>Q. (3.69) 

Similarly as in the regular case, we can choose the lens in such a way that /(S^nZ//) C T. 
Then P satisfies the jump condition ()3.36p for A G II5 H 

In order to have a good matching between P and on dlA, we define the analytic 
pre-factor E by 



i5^(A;e) = -^p(-)(A)(| (-.-^/V(A))^. 



(3.70) 



Using the definition (13:32]) of P^°°\ one checks directly that E is analytic in U. For 
A G and e~^/^(7i(A; x, tm)) e~^''''5'2(A; tm) in small complex neighborhoods of X, T E 
M, we have the matching 



P(A)P(~)(A)-1 = P~(A) (/ + iQcT3(-/(A))-i/2eV7 

(3 

as e — >• 0, with 

U = C/(e-6/7^i(A; X, t^), e-^/^92(A; i^)), (3.72) 

Q = Q(e-6/V(A; x, t„), e-^/'52(A; t™)), (3.73) 

= PF(e-6/V(A;x,t„),e-4/^52(A;t„.)), (3.74) 

y = y(e-6/V(A; X, t^), e-^/^52(A; *„)). (3.75) 

In the double scaling limit given by (jl.32p . we have by (|3.66p - p.67p that 
e-^l' gM-.^.i^) ^ e-6/^52(u'=;tm) ^ X, 

and from the definitions of g\ and (72, it follows that e^^/'^(7i(A; x, tm), e^^/'^g2(A; tm) 
lie in a small neighborhood of X^T if A is sufficiently close to u'^ . Summarizing, the 
matching (j3.7ip holds in the double scaling limit if we have chosen lA small enough. 
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3.6 Final transformation 5 i? 

We define 



R{X) 



fs(A)P(°°)(A)-i, asAeC\Z^, 
\s(A)P(A)-\ asXeU. 



(3.76) 



The outside parametrix has been constructed in such a way that the jumps of R are 
close to / as e — )• 0. On (0, +oo), the jump for S cancels out exactly with the jump of 
P(°°). On the other parts of S5 \ U, we have 



= l + 0{e), ase^O. 
For X gU n S5, we have 

RZHX)R+{X) = P^{X)vs{X)vp\X)PlH^)- 



(3.77) 
(3.78) 



(3.79) 



On one hand it follows from the construction of the parametrix that -P-(A) is uniformly 
bounded for A € n S^. On the other hand 



VSVp 



1 i{K — l)e e 
.0 1 



1 

l)e- 

K 

- \rr)e E 



K* + (l 



on Si n U, 
on E2 n u, 
on (n^O) nU. 



(3.80) 



Except for the 21-entry on (u^, 0), the exponentials in the above matrices are uniformly 
bounded on the jump contours inside U. The asymptotics (|3.22p and (j3.23p for k and 
r ensure that vsVp^ = I + 0{e). We obtain the following RH problem for R. 

RH problem for R 

(a) R is analytic in C \ S/j, with Sjj = (Sg U dU) \ (0, +00) as shown in Figure HI 

(b) R+{X) = R-{X)vfi{X), where the jump matrix vr is given by 



vr{X) 



'p(A)p(°°)(A)-\ asXGdU, 
l + 0{e), as X£T.R\dU ase^O. 



(3.81) 



(c) R{X) = 1 + 0{X'^) as A ^ 00. 

From the matching (j3.48p in the regular case where tm <t'^^ it follows that 

vr{X) = I + 0{e), ase^O, (3.82) 

uniformly for A G S/j, since is bounded on dlA. In the singular case, we use (|3.7ip 
to conclude that 



VR{X) = I + 0{e^/') 



(3.83) 
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Figure 4: The contour S/j after the final transformation S R. 



uniformly for A G S^, in the double scaling limit where e — t- 0, x — t- x'^, and — >• 
in such a way that (jl.32p holds. 

Following the general theory for small-norm RH problems [9], one shows that 
R{X) = I+— R_(s){vr{s)-I)^, forAGC\i?, (3.84) 
and that 

\\R--I\\LHJ:^) = 0{e^), (3.85) 

where 7 = 1 in the regular case as e — )• 0, and 7 = 1/7 in the double scaling limit in 
the singular case. 

4 Proof of the results 
4.1 Proof of Theorem [TH 



Note that (f3^ imphes that 

R: 



ii(A) = / + :^ + 0(A~2), as A ^00, (4.1) 



with 



R^ = R^{x,t^,e) = -— I R-{s){v{s) - I)ds. (4.2) 



R 



This gives 



Ri = -^ ( {R-{s) - I){v{s) - I)ds - / [vr{s) - I)ds = 0(6), (4.3) 

as e — )• by the Cauchy-Schwarz inequality. Differentiating the integral equation (j3.84p 
and (j4.3p in x, one similarly obtains the estimates 



= 0(e), = 0(e), (4.4) 

using the fact that ||w_R,a:||L2(j]^) = 0{e). Now by (j3.76p . we have that S{\) = 
i?(A)p(°°)(A) for large A, which implies by (15:351) that 

5ii(A) = l-i ^^'f' +0(A-^), as A ^00. (4.5) 
(-A)V^ 

By (j3.27p . we obtain 

u{x, tm, e) = u{x, tm) - 2edxRi,i2{x, tm, e) = u{x, tm) + O(e^), (4.6) 

which proves Theorem ll.2l in the case where x < xm + Cmtm- 
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Remark 4.1 If x = xm + Cmtrm the RH analysis remains the same, except for the fact 
that u = —1 and F'{u;x,tm) = +00. This imphes that the jump matrix for S decays to 
/ also near u, and there is no need to construct a local parametrix. If x > xm + Cmtm-, 
we have that —x + CmU™'~^tm + fR{u) = 0, with /r the inverse of the increasing part 
of no, but F{u]x,tjn) = —x + CmU'^~^tm + Il^u) does not vanish (in general). Because 
of this, we need to modify the G-function: 



„ {ri-X)y/ri-u J_i {'q-\)^u-r] 



n ( \ \ Vu — \ 

vr 

with 

ffRW 

'^W = / V A - uo{x)dx, for -1 < A < 0. 

This gives the additional jump condition 

Gm,+ W - Gm,-iX) = -2iT(A), for A G (-l,n). 
Then cj) takes the form 

2 3 

(j){\;x,tm) = -iriX) - Vu- XF{u;x,tra) + -{u - A) 2 F'(u; x, t^) 

A 



-±{u-Xf^F"{u;x,t^)-^ [ F'"(e;x,U(e-A)tde 
15 15 



With those modifications, the RH analysis can be carried on similarly as before, see 
also [n]. 

4.2 Proof of Theorem 11.41 

Using (I3.7ip . we can expand the jump matrix vr in fractional powers of e in the double 
scaling limit, 

vr{X) = 1 + eV^Ad) (A) + e2/7A(2) (A) + e^AaCs) (a) + 0{e^l'), (4.7) 

with 

AW(A) = iQ ■ {-f{X)y^/^P^'^\X)a3P^'^\X)-\ (4.8) 



(4.9) 



A(^)(A) = -i(-/(A))-ipM(A) ^2) P^^Hxr\ 

A(')(A) = -^(-/(A))-3/2pM(A) p(-)(A)-\ (4.10) 
for A G dU, and 

A(i) (A) = A(2) (A) = A(3) (A) = 0, for A G S^j \ dU, (4.11) 

since the jump matrices are equal to / up to an error of 0(e) on the other parts of the 
contour. Note that the functions A^^^ and A*-^) are meromorphic functions in U with 
simple poles at u'^, and that A'^'^^ is meromorphic in U with a double pole at u'^. Observe 
that the x-derivative of vr is not close to / on dU as e — )• 0, since the x-derivatives of 
U and Q cause multiplication with e'^/'^ by (l3:72]l - ([37731) . 
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Substituting (j4.7p in (j3.84p yields an asymptotic expansion for the RH solution R 
of the form 

R^X) = 1 + eV7i?{l)(A) + £2/7^(2) (A) + £3/7^(3) (;^) + 0(g4/7)_ (4_i2) 

Combining (j4.7p with (j4.12p and the jump relation i?+(A) = i?_(A)fij(A) gives the 
following relations for A G 

R^^\\) = i?W(A) + aW(A), (4.13) 
(A) = (A) + R^l^ (A) A(i) (A) + A(2) (A) (4.14) 
^ i?L^)(A) + i?L^^(A)A(2)(A) +iiL^)(A)A(i)(A) + a(3)(A). (4.15) 

In addition we know that i?(A) — )• I as A — )• oo, and consequently R^^\X) — )• for 
j = 1,2,3. So we have additive jump relations and asymptotic conditions for R^^\ 
R^'^\ and R^'^\ and it is easily verified that those conditions determine R^^\ R^'^\ and 
i?^3) uniquely. We obtain 

f^Res(A(|);.^), asA.C\Z7 

\^Res(A(i);u^) - AW(A), as A G Z^, ^ ^ 

i?(2)(A) = 

\^Res(i?«A«+A(2);n^), as A G C \ Z7, 



A 



:^Res(i?«A(i) + A(2);u^) -i?(i)AW(A) - A(2)(A), as A G ^. 



After a straightforward calculation we find using ()4.8p . (14. 9|) . and (I3.32p that, for A G 
C\Z7, 

i?(^)(A) = -g/V)-^/^V^fn (4.18) 



\-u^\0 
where 

Q = Q(e-6/V(n^;x,t^),e-^/^52(n^;tm)), 
U = C/(e-6A5l(n'=;x,^^),e-4/752(n^;^n^)). 
For the matrix R^'^\ using (|4.15p and (I4.16p - (l4.17p . one obtains 



R^'\X) = ^^^Res(A(2)(A);n^)+^^^^Res(A«(A);n=)+-^ 

A — A — A — 

+ --^^-^Res((A - n'=)A(3)(A); u^), for A G C \ Z7. 
(A — M^)^ 

We are interested only in the 12-entry of i?3, so it is sufficient to compute the entries 
in the first line of R'^^\u'^) and R'^'^\u^): 

R'n\u') = 0, 

(e-^^{u')Qx + e''/^g'^{u')QT Qf"{u')\ 

-^12 ) = \ - 



V V7V) 4/'(n- 

„(2). i>!Lrf £!M!^rr 

4?(«^) = o. 
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Consequently, 

^S'W = irb^+(A^^- '"^-^CXS, (420) 

where 



Z = {Q^Qx + UxQ - UQx + {V- W)x) 

^-4/7 / ( c\ 

+ {Q'Qt + UtQ - UQt + {V- W)t) 



2 



^{Q-^ + UQ + 3{V-W)), (4.21) 



8/'(u=)5/ 
and 

Compatibility of the small e-expansion (j4.12p with the large A-expansion (j4.ip learns 
us that 

RiM^, t, e) = -e^l'Q fiu'r^l^ + ^I'Z + 0[e^l'), (4.22) 



in the double scaling limit. Taking x-derivatives in (14. 3p . one justifies that we can 
formally differentiate ()4.22p to obtain 

-2e9.i?i,i2(x,t,e) = l^l'^ff^^U -2e'l' d.g^(u^)Zx 

2/'(n-)2 

(4.23) 

since Qx = U. Now we obtain by ([MSIl-dSSS]), the first equality in (gS]), and K2T\i 
that 

262^^/^ 



n(x,tm,e) = lie + l^-p- ) 



+ 



2(8A;)6/7 

52 (ii^) /2(m-l) 3/" 



(8fc)4/7 V 4 / 



7) {2UxQt + ^UUt + \uxxt) + 0(e5/^). 



In the derivation of this expansion we have used (j3.54p - (j3.56p . For the third term, we 
also used the identity (jl.lOp . The term in the third line is also of order 0(e^/^) because 
g2iu^) is of order 0{e^^'^). The expansion for u{x,tm,^) can be written in terms of 
U exclusively by substituting the formula (|3.54p for Q and Qt = — C/^/2 — Uxx/'^2. 
Substituting the values for f"{u^), gi{u'^), and 52 ("u"^) given in (|3.65p . (j3.66p . and (|3.67p . 
we find (jl.25p . with the constants ci, C2, C3 given by ()1.28p - (jl.30p . This proves Theorem 

m 
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